Astronomy & Astrophysics manuscript no. magneticspots 


©ESO 2011 


August 11,2011 





Magnetic spots on hot massive stars 



M. Cantiello & J. Braithwaite 



Argelander-Institut fur Astronomie der Universitat Bonn, Auf dem Hiigel 71, D-53121 Bonn, Germany 
e-mail: cantiello@astro . uni -bonn . de 



ABSTRACT 



o 

(N 
< 



in 

6 

-I— > 



> 
o 
m 
O 

(N 
00 

o 



X 



Context. Hot luminous stars show a variety of phenomena in their photospheres and winds which still lack clear physical explanation. 
Among these phenomena are photospheric turbulence, line profile variability (LPV), non-thermal emission, non-radial pulsations, 
discrete absorption components (DACs) and wind clumping. Cantiello et al. (2009) argued that a convection zone close to the stellar 
surface could be responsible for some of these phenomena. This convective zone is caused by a peak in the opacity associated with 
iron-group elements and is referred to as the "iron convection zone" (FeCZ). 

Aims. Assuming dynamo action producing magnetic fields at equipartition in the FeCZ, we investigate the occurrence of subsurface 
magnetism in OB stars. Then we study the surface emergence of these magnetic fields and discuss possible observational signatures 
of magnetic spots. 

Methods. Simple estimates are made using the subsurface properties of massive stars, as calculated in ID stellar evolution models. 
Results. We find that magnetic fields of sufficient amplitude to affect the wind could emerge at the surface via magnetic buoyancy. 
While at this stage it is difficult to predict the geometry of these features, we show that magnetic spots of size comparable to the local 
pressure scale height can manifest themselves as hot, bright spots. 

Conclusions. Localized magnetic fields could be widespread in those early type stars that have subsurface convection. This type of 
surface magnetism could be responsible for photometric variability and play a role in X-ray emission and wind clumping. 



1. Introduction 

There is growing direct and indirect evidence that magnetic 
fields are present at the surface of massive stars. Direct measure- 
ments have been obtai ned using t he Zeeman effect for O and 
early B stars (see e.g. lPetitll2010l for a review). As the num- 
ber of known magnetic massive stars is small, these observa- 
tions have not yet been able to answer the question of whether 
magnetism in these stars is fundamentally different from that 
in the intermediate-mass stars, i.e. late B and A stars (the so- 
called Ap and Bp stars; see e.g. iMathvsl J20 09) for a review, or 
Mestel & Landstreet] d2005l) for a broader review of stellar mag- 
netism), or just a continuation of the same phenomenon. Ap and 
Bp stars - accounting for some percentage (~ 5 - 10%) of main- 
sequence A and late B stars - display a magnetic field which 
is static, in the rotating frame (these stars appear to rotate as 
solid bodies), large-scale, and with strengths of 200G to lOOkG. 
These fields are thought to be 'fossil' remnants left over from 
formation: magnetohy drostatic equilibria which survive for the 
entire stellar lifetime ([Cowling! 1 1945b Braithwaite & Nordlundl 
2006; iDuez et al.ll2010l) . However, since direct detection meth- 
ods are biased toward the discovery of large scale magnetic 
fields, small-scale fields generated by a contemporary dynamo 
may be present in the rest of the population. Whilst in A stars 
even smal l-scale fields w ould have to be rather weak to escape 
detection (Petit et al.ll2010h . the detection limits on small-scale 
fields in more massive stars are much higher (e.g. lSchnerr et alj 
l20Q8h . 

In addition there is a variety of indirect evidence for such 
magnetic fields at the surface of OB stars, coming from ob- 
servations of surface and wind phenomena. For example, dis- 



crete absorption components (DACs) and line profile variabil- 
ity (LPV) have been often associated with the presence of sur- 
face magnetic fi elds, even if non radial pulsation (NRP) is also a 
possibility (e.g., Cranmer & Owockill996l ; lFullerton et ai7 l996: 



al 1997; 



Henrichs et al ■1120051) . Interestingly these phe- 
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nomena seem to be ubiquitous among early type stars (e.g. , 
[ Howarth & Prinidl989tlKaper et alJll996l;lFullerton et al.|[ l996: 
iPrinia et aUl2002l) . Therefore one may wonder if magnetic fields 
in massive stars are more widespread than previously thought. It 
is important to stress that small scale magnetic fields of ampli- 
tude ~10...10Q G are in principl e able to affect the wind of OB 
stars (ud-Doul a & Owockll2002h . Moreover, theoretical consid- 
erations point to a difference between massive and intermediate- 
mass stars. In the latter, a dynamo would have to operate either 
in the radiative envelope, powered by differential rotation ( Spruit 
2002) which itself would have to be sustained by some physical 
process; or in the convective core, in which case the resulting 
field would h ave to make its way upwards through the radia- 
tive envelope dMossll989HMacGregor & Cassinellil2003l) which 
seems to take place only on an embarrassingly lo ng timescale. In 
contrast to this, recently ICantiello et alJ d2009l hereafter C09) 
studied the properties of convective regions in the envelopes 
of hot, m assive stars, cause d by an opacity peak associated 
with iron (Igles ias et al ] l!992h . These convective zones (here- 
after FeCZ), which are not present in intermediate-mass stars, 
could host a dynamo, producing magnetic fields visible at the 
surface. 

The reality and importance of subsurface convective lay- 
ers in early-type stars has been recently stressed by observa- 
tions of stellar surface properties. The idea that the observed mi- 
croturbulence is caused by the pr esence of the FeC Z (C09) is 
supported by recent observations (Fra ser et al.ll2.Q10L Przybilla, 
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priv.comm.). In the S un stochastic convective motions excit e 
non-radial pulsations (lUlrichl Il970t iLeibacher & Steinl Il97ll) . 
In analogy with the Sun, it has been suggested that solar- 
like o scillations are excited in hot, massive stars by the FeCZ 
(C09, iBelkacem et alllloToh . Intriguingly, non-radial pulsations 
that are compatible with stochastic excitation from subsurface 
convection have recen tly been found by COR OT in both a B 
(Belk acem et al.ll2009l) and an O star dDegroote et al.ll2010h . but 
see also lBalona et al.l (1201 ll) . 

C09 pointed out that the FeCZ could host dynamo action 
very close to the surface. Simulations of turbulent convection in 
the presence of shear and rotation show that dynamo action is 
possible and can produce magnetic fields reaching equipartition 



on scales larger than that of the convection (e.g ., Kapvla et"al] 



2008; Cantiell o et all 120101: iGuerrero & Kanvlall201 lh . In this 
paper we consider the emergence of magnetic fields produced in 
a subsurface convective layer. In the following section we briefly 
review the properties of the FeCZ, before examining in Sect. [3] 
the magnetic fields produced within them and mechanisms for 
those fields to reach the surface. In Sect. [4] we look at observ- 
able effects on the surface, and finally we discuss the results and 
conclude in Sect.[5]and[6] 



2. Subsurface convection 

The occurrence and properties of the iron convection zone were 
studied in detail by C09. Here we summarize their findings and 
describe some of the fundamental properties of the outer regions 
of hot massive stars. The evolutionary calculations used in this 
paper are the sam e as in C09. We refer to that paper and to 
Brott et al] d201 lh for a detailed discussion of the models and 
the stellar evolution code with which these are computed. 

The peak in the opacity that causes the FeCZ occurs at 
around 1.5 xlO 5 K. In the envelope of a massive star this region 
remains relatively close to the surface during the whole main 
sequence. When the star expands to become a supergiant, the 
surface temperature decreases and the convective zone moves 
downwards in radius, further from the surface. For photospheric 
temperatures < 10 4 K, partial ionization of H and He drive con- 
vection directly at the surface. Therefore we limit the discussion 
of the properties of the FeCZ to spectral classes O and B (and 
partially A), which lack this additional convective layefl 

The presence of subsurface convection depends also on the 
luminosity and metallicity of the star. At solar metallicity the 
FeCZ is present in models above ~ 10 3 2 L Q , while at the metal - 
licities of the LMC and SMC the thresholds are ~ 1O 3 9 L and 
~ 10 4 2 L respectively. Atomic diffusi on and radiative acc elera- 
tion can in principle lower these limits (Richer et al. 2000). 

Due to the very low densities in the outer layers of early 
type stars (p ~ 10~ 8 g cm -3 ), these convective regions contain 
very little mass despite their large radial extent. In Table [TJ we 
show the properties in the outer layers of Galactic (solar metal- 
licity) O/B stars above the luminosity threshold. Values for the 
size of the convection zone and of the radiative layer above are 
shown for models of a 20 and a 60 M Q at ages 6.41xl0 6 yr 
and 2.37xl0 6 yr respectively. These correspond to the calcula- 
tions discussed in detail in C09, where is possible to see how the 
extension and location of subsurface convection changes during 
the MS evolution (See their Fig. 2 and 3). 

In the FeCZ the transport of energy by convective mo- 
tions is relatively inefficient: radiation dominates and transports 



most of the total flux (e.g. about 95% in the 60 M model of 
Table [T]). This is because the density is very low and the mean 
free path of photons correspondingly long. In this situation con- 
vection is significantly superadiabatic, and the actual gradient 
V = d In T/d In P has to be explicitly calculate d from the mixing 
length equations (e.g. Kippenhahn & Wei gerH fl 9901) . C09 found 
velocities ranging from 1 to 10's of kms" 1 in the FeCZ of OB 
stars, where the sound speed is usually of the order lOOkms -1 . 

In addition, it is informative to estimate the values of the 
diffusivities in the FeCZ. As usual in non-degenerate stars, the 
thermal diffusivity is the largest: 
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where a, c p and k are the radiation constant, the specific heat 
at constant pressure and the Rosseland mean opacity. In the 20 
and a 60 M models^ is around 5 x 10 17 and 8 x 10 19 cm 2 s _1 
respectively. The momentum diffusivity (kinematic viscosity) is 
(in contrast to lower mass stars) also dominated by photons, and 
is given in this case by 
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where in the 20 and a 60 M models it is approximately 6 x 10 9 
and 10 12 cm 2 s _1 . Fi nally, the mag netic diffusivity as given by 



Spitzer's expression (Spitzer 1962) 
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(3) 



1 The discussion in Sect.f3]is more general, and can be applied to any 
radiative star with a convective region below the surface. 



which is roughly 10 cm s in both models (assuming a fully 
ionized H plasma and a value In A = 10 for the Coulomb log- 
arithm). Note that photons have a negligible effect on the mag- 
netic diffusivity. The magnetic Prandtl number (Pr m = v/rj) is 
consequently 10 5 to 10 7 . We can compare these numbers to the 
values usually discussed in the solar dynamo (Pr m ~ 10~ 2 in the 
tachocline) and in the galactic dynamo contexts (the ISM has 
typically Pr m ~ 10 12 ). 



3. Magnetic spots 

The occurrence of convection zones close to the surface of hot 
massive stars opens a new scenario. Magnetic fields could be 
readily produced by dynamo action and reach the surface via 
magnetic buoyancy (C09). Below we discuss this hypothesis. 



3. 1 . Dynamo Action 

In an astrophysical plasma a dynamo is a configuration of the 
flow which is able to generate a magnetic field and sustain it 
against ohmic dissipation. Depending on the scale of the result- 
ing magnetic field respect to the scale of kinetic energy injec- 
tion, dynamos are usually divided into small and large scale. 
In large scale dynamos the field has correlation length bigger 
than the forcing scale in the flow, while small scale dynamos re- 
sult in magnetic fields with correlation scale of order or smaller 
than the forcing scale. Small scale dynamos can occur in non- 
helical turbulent flows, while anisotropic flows (e.g. shear flows) 
are required for a large scale dynamo. In the mean field ap- 
proach, large scale dynamos are often divided into o-Q and a 2 , 
depending on the role of the D. and a effect in regenerating 
the toroidal and poloidal components of the field. We refer to 
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Table 1. Properties of the outer layers in a 20 M and 60 M model main-sequence stars of solar metallicity at ages 6.41xl0 6 yr 
and 2.37xl0 6 yr respectively, corresponding to about 90% and 80% of the main-sequence lifetime. 
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" Radial coordinate of the top of the FeCZ. 

Radial extension of the FeCZ. 
c Pressure scale height at top/bottom of the FeCZ. 
d Maximum of the convective velocity inside the FeCZ. 
' Density at v c . 

Mass contained in the convective region. 
s Mass in the radiative layer between the stellar surface and the upper boundary of the convective zone. 

Convective turnover time, T lurn := Hp/v c . 
' Time that a piece of stellar material spends inside a convective region, T conv := AM^z/M. 



Brandenburg & Su bramanianl d2005l) for a very nice review of 
astrophysical dynamos. 

The FeCZ is a turbulent layer close to the surface of a mas- 
sive star and as s uch could be the site of a small scale dynamo 
(e.g.. lMosslll994l) . Moreover main sequence massive stars usu- 
ally rotate rapidly. A typical equatorial rotationa l velocity of 
150 k ms' 1 (typical for Galactic B-type stars, e.g. iDufton et al.l 
2006) corresponds to a rotational period of the order of the 
convective turnover timescale (Rossby number is in the range 
1 ... 10), therefore an efficient a 2 or aQ-dynamo could be possi- 
ble. Assuming magnetic fields at equipartition with the kinetic 
energy of the convective motion, gives magnetic fields up to 
~ 2kG. This is supported by simulations of turbulent convection 
in the presence of rotation and shear, which show dynamo excita- 
tion with magnetic fields rea ching equipartition on scales larger 
than t he scale of convection (Kap vla et alJl2008t ICantiello et al.1 
l201Ch . 

Hence dynamo action in the FeCZ could depend on parame- 
ters like the stellar rotation (as measured by the Rossby num- 
ber) and the shear profile in the region of interest. The scale 
of the magnetic field depending on the type of dynamo occur- 
ring in the relevant layers. The dynamo may also be affected by 
a (presumably fossil) large-scale magnetic field penetrating up- 
wards from the radiative zone below (see Sect. |5). However it 
is interesting to note that many of the photospheric and wind 
phenomena observed in hot massive stars are ubiquitous (e.g. , 
Howarth & Priniall989tlkaper et al.lll996UFullerton et al.ll996t 
Prinia et al.l 20021 in contrast to fossil fields which have been 



B(G) 



found at the surfaces of some (as yet badly measured) fraction 
of massive stars. As it is difficult to predict the exact rotational 
properties of the plasma inside the FeCZ, here we will not focus 
on a detailed study of subsurface dynamo action. As a prelimi- 
nary stud y of subsurface mag netism we will just follow the re- 
sults of Cantiello et al. (2009, 2010), and assume that magnetic 
fields at equipartition are produced in the FeCZ. In Fig. [TJ we 
show the expected maximum magnetic field inside the FeCZ, 
assuming equipartition of magnetic and kinetic energy: 



B eq 1 2 

to = 2 pV » 



(4) 



where we adopted the maximum value of the pv\ inside the FeCZ 
(as computed in the non-rotating models of C09). Due to the 
higher power dependency on the velocity, the location of the 
maximum of pv\ always roughly corresponds to the maximum 
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Fig. 1. Values of maximum magnetic fields (in gauss) in the 
FeCZ, as a function of the location in the HR diagram. This 
plot is based on evolutionary models between 5 M and 120 M 
for solar metallicity (some evolutionary tracks shown as dotted 
lines). The amplitude of magnetic fields is calculated assuming 
equipartition with the kinetic energy of convective motion. The 
full drawn black line roughly corresponds to the zero age main 
sequence. The two star symbols correspond to the location of the 
20 and 60 M models discussed in the text and in C09. 



of v c . Values of magnetic fields calculated using the average con- 
vective velocity and density typically differ by less than 30%. 

3.2. Magnetic buoyancy 

There are various possible mechanisms which can bring a mag- 
netic field generated in the FeCZ through the overlying radiative 
layer to the surface, and it is a useful exercise to compare their 
timescales. As the title of this section suggests, it turns out that 
among the mechanisms studied here, the dominant one is mag- 
netic buoyancy. 

The magnitude of the mass loss from massive stars is such 
that material resides in the outer, radiative layer for just a short 
time. Assuming a rate of ~ 1O~ 7 M yr -1 , this time is ~ 1 yr. 
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Magnetic field could also be brought to the surface by Ohmic 
diffusion. The timescale over which it produces effects over the 
relevant radial length scale, i.e. the pressure scale height Hp, is 



TOht. 



Hi i H P \ 2 ( T \ 3/2 v 



where tj is the magnetic diffusivity. We shall now see that this 
timescale and that of advective transport are much longer than 
the buoyancy timescale. 

Buoyancy causes magnetic features to rise through a radia- 
tive zone because magnetic field provides pressure without con- 
tributing to density. Since the sound and Alfven timescales are 
shorter than the thermal/buoyancy timescale (the buoyant rise is 
subsonic and sub-Alfvenic), we can assume that a magnetic fea- 
ture is in pressure equilibrium with its surroundings. Calling the 
sum of gas and radiation pressure inside and outside the feature 
Pi and P e we have 

P e = Pi + Praag, (6) 

where the value of the magnetic pressure P mag depends on the 
geometry of the feature. For instance, in a self-contained struc- 
ture it is B 2 /24n, and in a flux tube of fixed length it is B ax /87r 
where B ax is the component along the axis of the tubeQ. Note 
the implicit assumption that the size of the structure / is much 
smaller than the scale height Hp; it is quite likely (see below) 
that in fact / » Hp but this will make a difference just of fac- 
tors of order unity. Now, in the absence of thermal diffusion 
the feature reaches an equilibrium where p\ = p e , made pos- 
sible by the internal temperature being lower than the external 
T\ < T e . With the addition of a small thermal diffusion, the fea- 
ture absorbs heat from its surroundings and rises quasistatically 
upwards through surroundings of increasing entropy. The speed 
of this rise is given by: 



2xHp 
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(7) 



where B = P e /P mas , I is the size of the feature, V and V a d have 
their usual definitions and take values V = V la d ~ 0.24 and 
V a d ~ 0.25 in the radiative zones between the FeCZ and the 
photosphere, and^- is the thermal diffusivity given in (fjji. 

If on the other hand thermal diffusivity is large, the feature 
rises so fast that the speed is limited by aerodynamic drag. In this 
regime the speed is independent of thermal conductivity, and we 
can make the approximation that 7j = T e ; note that unlike above, 
in this regime we have Ap = p e - p; + 0. The buoyant force as 
given by Archimedes' principle is balanced by the drag force, 



1 



:C d Ap e v 2 d = VgAp 



(8) 



where A and V are the cross-sectional area (projected onto a hor- 
izontal plane) and volume of the magnetic feature, and Cd is the 
drag coefficient whose value depends on geometry. This can be 
rearranged to 

I 2V \P 



^[CaHpAj p, 



mag 



(9) 



Since the length scale of the feature I « V/A, the term in brackets 
is approximately equal to ////p0 The remaining part is approx- 
imately the square of the Alfven speed, which one can alterna- 
tively express in terms of the sound speed c s and the ratio of total 



2 See Braithwaite (2010) for a discussion of this point and a review 
of the stability of flux tubes. 

3 If the feature is spherical, d ~ 0.75 tCh urazov et al.ll2~00lh . 



to magnetic pressure, B: 
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where ranges in the numbers come from the difference between 
20 and 60 M models and from differences between top and bot- 
tom of the surface radiative layer. It is clear that unless I » Hp, 
for which the approach above breaks down anyway, or the field 
is implausibly weak (and therefore undetectable and unimpor- 
tant), Utherm is larger and the rise is therefore limited by drag. 
For a magnetic feature with / ~ Hp and B ~ 100 (a conserva- 
tive estimate), the time taken to rise one scale height is of order 
2 hours in both models, much shorter than the advective and 
Ohmic timescales. 

This is in contras t to the thermal-diffusion-limite d regime 
which was studied by iMacGregor & Cassinellil (120031) to check 
the possibility that magnetic fields generated in a convective core 
could rise to the surface. In that case the timescale appears to ex- 
ceed the main seque nce lifetime unless very smal l-scale features 
are used. Moreover Mul lan & MacDonaldl (12005) argue that the 
inclusion of compositional gradients may suppress the buoy- 
ancy, making this scenario unlikel)Q. 

The situation in massive star surface layers is in some sense 
more similar to that of flux tubes in the solar convection zone 
in that the rise is limited by drag forces, although the drag force 
in the solar CZ does not scale simply as v 2 but depends also on 
the convective motion. A stationary flux tube in the solar CZ 
experiences a net downwards aerodynamic force once the av- 
erage of upwards- and downwards-moving convective cells has 
been taken, with the result that there is a field-strength thresh- 
old below which the tube moves downwards despite its intrinsic 
buoyancy. 

Once a magnetic field element reaches the surface of the star 
it would be useful to estimate its lifetime f spot . This is not trivial, 
as the photosphere of OB stars is complicated by the presence of 
strong winds and, likely, turbulence and shear. However we can 
put some limits on the lifetime. 

As we discussed at the beginning of this section, if the mate- 
rial in the radiative layer above the FeCZ has a mass AM, it will 
be removed by mass loss in a time AM/M. This sets an upper 
limit for the time a magnetic element, which is not anchored to 
the convective zone, can exist at the surface: this time is about 
8 years for the 20 M and 3 months for the 60 M model. Of 
course this holds true only as long as the field is not annihilated 
by dissipative processes like magnetic reconnection. If the mag- 
netic element is anchored to the convective zone, it can exist at 
the surface as long as the stellar wind does not remove all the 
mass contained in both the radiative layer and the FeCZ. It turns 
out that this limits f spot to be less than 50 years (20 M Q model) 
and 4 years (60 M model). We refer to Table 1 in C09 for the 
data used for these simple estimates. These are solid upper lim- 
its, and it is likely that such features would survive for a much 



4 Mullan & MacDonald ( 2005) also discuss the possibility of buoy- 
ant rise of magnetic fields generated in the radiative zone by the Tayler- 
Spruit dyn amo. However it is not clear how thi s dynamo process actu- 
ally works (Braithwaite 2006; Zahn et al. 2007). 
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Fig. 2. Schematic representation of the effect of an emerging 
magnetic element at the surface of a hot massive star. A magnetic 
field (dashed line) rising from the subsurface convection zone 
threads the radiative layer and reaches the stellar photosphere. 
The solid, grey line shows the location of the stellar surface, as 
defined by the value 2/3 for the optical depth. Notice that inside 
the magnetic field this line reaches deeper (and hotter) layers of 
the star, as compared to regions not affected by the field. 

shorter time, owing to turbulence, shear and magnetic reconnec- 
tion. As a lower limit for f spot we could assume the time a feature 
of size I takes to cross the photosphere while rising with a ve- 
locity Vd ra g. For I ~ Hp this gives a timescale of the order of 
hours. 

Finally, note that it is not inconceivable that the magnetic dy- 
namo does not saturate properly because of the finite time the gas 
spends in the FeCZ. However, we consider this unlikely, since 
the convective turnover times (~ hours - a day) appear to be very 
much shorter than the aforementioned mass-flux timescales of 
50 and 4 years. 

4. Observable effects 

Let us first estimate the likely field strength at the photosphere. 
Whilst one might expect it to be proportional to the field strength 
produced in the convective layer, the difference in field strength 
of a magnetic feature between the top and bottom of the radiative 
layer depends on its geometry. In a self-contained magnetic fea- 
ture (a 'blob' or 'plasmoid') the field strength scales as B oc p 2 ^ 3 , 
while a horizontal flux tube of fixed length has B ax oc p and 
Bh oc p 1 ^ 2 where B ax and B^ are the axial and 'hoop' compo- 
nents. Knowing that P ~ p 4 ^ 3 (approximately, since V«l /4) in 
the radiative layer, we see that these two scenarios give B - const 
and B oc p -1 / 2 respectively. In constrast, if the central section of 
a flux tube rises to the surface while its ends are still in the con- 
vective layer, plasma can flow downwards along the tube and 
the field strength at the surface can in principle reach equipar- 
tition with the surrounding (gas plus radiation) pressure, as we 
see happening in the solar CZ. Assuming equipartition magnetic 
field in the convective layer gives up to « 2 kG at the base of 
the radiative layer, and B as 100. The density contrast across the 
radiative layer is around 50 and 100 in the 2OM and 6OM mod- 
els respectively, which gives field strengths at the photosphere of 
40- 150 G and 20- 100 G respectively in the blob and horizontal 
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Fig. 3. Minimum values of expected surface magnetic fields (in 
gauss) as a function of location in the HR diagram. This plot 
is based on evolutionary models between 5M and 120 M Q 
for Galactic metallicity (some shown as dotted lines). Surface 
magnetic fields are calculated scaling the equipartition fields 
in the FeCZ (see Fig. Q} with p 2 ^ 3 (see discussion in Sect. |4). 
The full drawn black line roughly corresponds to the zero age 
main sequence. Cir cles correspond to the observed stars in the 
ICohen et alJ (1 19971) sample, with symbol size coding different 
range of L x /L\, Q i, as explained in the legend. The two star sym- 
bols correspond to the location of the 20 and 60 M models dis- 
cussed in the text and in C09. 



tube cases. The photospheric field could be significantly greater 
if a tube arches upwards and fluid is allowed to flow along it - 
the limit should then be of the order of equipartition with the 
photospheric pressure (about 300 G) or somewhat more, since 
the photosphere (i.e. the r = 2/3 level) in the tube is lower. 

In Fig. [3] we show the expected amplitude of surface mag- 
netic fields emerging from the FeCZ, assuming a scaling B oc 
p 2 ^ 3 . We argue that this represents a lower limit for the mag- 
netic field that should be found at the surface of hot massive 
stars, if equipartition of kinetic and magnetic energy is reached 
in their FeCZ. It might seem plausible that magnetic features ris- 
ing through the radiative layer partially annihilate each other on 
the way up. However, if these magnetic bubbles have a size at 
the bottom of the radiative layer comparable to the local scale 
height, they only have to rise through a distance about three 
times their own size to reach the photosphere, expanding as they 
do so. They rise at approximately the Alfven speed, and since 
any reconnection would presumably take place at some fraction 
of the Alfven speed, there is limited time available to destroy 
much of the original flux. It seems more likely that significant 
reconnection does take place above the photosphere. 

While these fields appear to have very small amplitude, their 
role at the stellar surface should not be underestimated. In Fig. [3] 
the ratio of the magnetic t o wind pressure (rj tss B^Jt^/Mv^, 
lud-Doula & Owockl |2002|) is of the order of unity over the 
whole region where magnetic fields are present. This means that 
already such low magnetic fields can alter the dynamic of the 
stellar wind, for example contributing to seed instabilities in the 
flow. 
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4.1. Surface bright spots and photometric variability 

Once a magnetic field element reaches the stellar photosphere, 
one would like to know what effect this can have on observable 
properties of the star. In the Sun, magnetic fields of amplitude 
up to ~ 3 kG emerge at the surface as sunspots. The main effect 
of such magnetic fields is to locally inhibit convection. As the 
energy flux below the surface of the Sun is mainly transported 
by convection, this produces a decrease of the flux at the sur- 
face, resulting in a visible dark spot. In contrast, the hot massive 
stars considered here have radiative surfaces and the effect of an 
emerging magnetic field is very different. 

Under the simple assumption of magnetohydrostatic equi- 
librium we can make use of Eq. [6] Assuming again that also 
the gas temperature inside and outside is the same (T e = T,), 
we obtain once more that the density inside the magnetic ele- 
ment is lower than outside. Since we are at the stellar surface, 
this makes the element more "transparent", i.e. photons can es- 
cape from deeper down the star compared to outside the mag- 
netic flux. This is somehow similar to what occurs in the solar 
faculae, which are regions of emerging magnetic fields with a 
scale much smalle r than the scale of convection. In that case 
iKeller et all (12004 found, using MHD calculations, that "The 
opacity in the magnetic flux concentration is strongly reduced 
owing to its low density and temperature and thus provides a 
clear sight straight through the magnetic field onto the adjacent 
nonmagnetic granule". Their calculation s basically verified t he 
explanation for solar faculae proposed by Spruit (II 97611 19771) . 

Photons emerging from deeper down the star are emitted 
from regions where the temperature is higher. The tempera- 
ture stratification, assuming that the magnetic pressure is much 
smaller than the gas + radiation pressure (i.e. /3 » 1), is given 
by the radiative gradient V ra( j 



V rad = 



dlnT 
dlnf 



Therefore it is easy to estimate the temperature difference be- 
tween the region of the star visible through the magnetic element 
and outside of it: 



AT 
AP 



TAQnT) J v 
PA(lnP) ~ P rad 



Ar 



Vrad 
P 



(13) 



This means that regions with an emerging magnetic field will 
look hotter. In the case of the FeCZ, assuming the minimum 
value of magnetic fields shown in Fig. [3] preliminary estimates 
give temperature differences up to a few hundred K. In the case 
of a surface magnetic field in equipartition with the photospheric 
pressure (gas plus radiation) equation[T3]is not strictly correct, as 
the magnetic field can affect the radiative gradient. Nevertheless, 
assuming for a moment that the correction to this relation is 
small, one expects temperature differences up to a few kK. 

Since the stellar luminosity is related to the effective temper- 
ature through the relation L = 4-7tR 2 ctT 4 ', magnetic fields in stars 
having radiative surfaces will produce bright spots. We can esti- 
mate the local luminosity contrast between the magnetic bright 
spot and the non magnetic regions of the stellar surface 



ALi r 



(AT\ 4V ra 



P 



(14) 



Temporal changes in the observed stellar luminosity might be 
expected, as a result of the appearance/disappearance of such a 



bright spot. In a rotating star the luminosity would also decrease 
(increase) when the spot moves out of (into) the visible part of 
the stellar disc. It is difficult at this stage to give a firm estimate 
of the total expected change in luminosity, since this will depend 
not only on the intensity of the magnetic field, but also on the 
geometry and number of magnetic spots. Also, it is not clear ex- 
actly where the extra luminosity in the spots should come from, 
in terms of underluminous areas. Nevertheless, if we assume a 
circular spot of radius r, with filling factor / = (r/R) 2 , where R 
is the stellar radius, then the effect on the stellar luminosity can 
be estimated if r is known. A simple ap proach, based on prel im- 
inary models of subsurface convection (Cantiello et al1 l2010l) . is 
to equate r to the pressure scale height in the FeCZ. These simu- 
lations show dynamo action producing magnetic fields on scales 
comparable to or larger than the scale of convection, i.e. the local 
pressure scale height Hp. We obtain: 



AL 



4V r 



rad 



P 



(15) 



which, for typical values of subsurface convection, gives AL/L ~ 
10~ 5 . This is the relative change in luminosity caused by the ap- 
pearance/disappearance of one magnetic spot of about 100 G 
with size 0.2 R Q at the surface of a hot massive star. Recall that 
these stars have radii around 10R Q . 

Of course spots with stronger magnetic fields and bigger fill- 
ing factors are certainly possible. Such spots would lead, most 
likely, to larger changes in the integrated light from the star. An 
estimate of the amplitude of the variability as function of filling 
factor and surface strength is much more difficult in this case, 
as one can no longer consider the magnetic element thin (in a 
thermal sense). MHD calculations of magnetic flux tubes in a 
radiative environment are required to study the details of this 
problem. 

Photometric variability at the micro level (10~ 3 magni- 
tudes) and with time scales of the order of days seems to be 
widespread in O stars (Balona 1992). U sing COR OT to perform 
high-precision photometry, Degr oote et alJ d2010h found solar- 
like oscillations in a young O-type star (HD46149). These are 
modes with a finite lifetime that are be lieved to be excited by 
the presence of subsurface convection (Belkac em et alJ l2010). 
Intriguingly, this system also shows a low-frequency variabil- 
ity with time scale of days and amplitude of order 10~ 3 magni- 
tudes. The authors could not relate this variability to any type 
of pulsations. Instead they argue that this could be the signa- 
ture of spots, stellar winds or chemical inhomoge n eities modu- 
lated b y the stellar rotation (Degroote et al.l 12010 ). Papics et al. 
d2oTTh obtained high-precision photometric observations of the 
B0.5IV star HD51756 with COROT. While no solar-like oscil- 
lations were identified in this case, they found cusp-like fea- 
tures in the light-curve. These features have amplitudes of or- 
der 10~ 3 magn itudes, and recur o n a time scale of days. Similar 
to HD46149, [Papics et al] (1201 ll) propose as a possible source 
for this variability photospheric features like spots or chemical 
inhomogeneities, modulated by stellar rotation. 

4.2. X rays 

OB s tars have been found to emit X rays (Harnde n et al.l 
119791). O stars follow th e relation L x /L ho i ~ 10" 7 (see 
e.g. Pallavicini et al ] 119811) while B star X -ray luminosities 
show a lot of scatter (IMeurs et all 1 19921: [Berghoefer et al. 
1 1997b . Among possible physical me chanisms for the produc- 
tion of X rays are magnetic coronae dCassinelli & Qlsonlll979l : 
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IWaldron & Cassinellil 12007]. [2 009) and the intrinsic instability 
of line-driven winds (iLucv & Whitdll980t lOwocki et alJll98a 
lFeldmeierfll995h . 

The L x /Lboi relation only applies for Lboi > 10 38 erg s _1 , i.e. 
earlier than B0-B1 type stars. It breaks down at lower luminosi- 
ties, with the r atio L x /Lboi two orders of magnitude s maller at 
B3 than at Bl dCassinelli et al]ll994tlCohen et : al.lll997l). This is 
shown in Fig. [3] where we plot the data from lCohen et al.l (11997) 
together with the theoretical expectation for the minimum val- 
ues of surface magnetic fields (assuming dynamo action in the 
FeCZ). While there could be different explanations for such a 
correlation, it is interesting to note that in the sample stars with 
Lx/Lboi > -7 correspond to regions of the HRD wi th larger val- 
ues of predicted surface magnetic field. Recently Gagn e et alj 
d20Tlh have studied with Chandra the X-ray emission from OB 
stars in the Carina Complex. They also found that most of their 
O stars follow the L x /L\, \ relation, while this breaks down for 
B stars. Interestingly they identify a group of B stars with high 
X-ray emission that cannot be explained by a distribution of or- 
dinary coronal, pre main-sequence companions. A possibility is 
that in O stars wind embedded shocks dominate and are respon- 
sible for the L x /Lboi relation, while the observed X-ray luminos- 
ity in B stars is mainly due to the presence of subsurface gener- 
ated magnetic fields. 

5. Discussion 

We showed that, among advection, magnetic diffusion and mag- 
netic buoyancy, the most likely process that can bring to sur- 
face magnetic fields generated by dynamo action in the FeCZ 
is magnetic buoyancy. However there could be other ways for 
a magnetic field to escape the subsur face convection region and 
reach t h e stellar surface. For exa mple JWarnecke & Bra ndenburg 
(2010); Warne cke et al.l d201 ll) studied the magnetic flux pro- 
duced by a turbulent dynamo in Cartesian and spherical geom- 
etry, respectively. They found that magnetic flux can rise above 
the turbulent region without the need of magnetic buoyancy. This 
appears to be related to the release of magnetic tension, which 
leads to the relaxation and emergence of the field. Convective 
overshoot might also contribute to the transport of magnetic flux 
out of the FeCZ. However, given the fact that the FeCZ usually 
sits a few pressure scale heights from the stellar surface, this ef- 
fect is likely marginal. 

Some fraction of massive stars display large-scale, appar- 
ently fossil fields of strengths (1 -3 kG) comparable to the FeCZ 
equipartition strengths derived above. A sufficiently strong field 
will suppress convection, forcing an increase in temperature gra- 
dient so that the entire energy flux can be transported radiatively. 
However, it is not obvious where the field strength threshold 
should be; one might naively expect convection to be suppressed 
by a field of gr eater energy density tha n the convective motion , 
but the work o f lGough & Tavlen 1 (fl966h : lMoss & Tavloil dl969h : 
lMestell(ll970l) suggests that in that case the temperature gradient 
would simply steepen further above the adiabatic gradient until 
convection resumes, and that to suppress convection completely, 
a field above (approximate) equipartition with the thermal en- 
ergy would be required. However, this work considers a situation 
where the energy flux is almost entirely convective, such as in the 
bulk of the solar convective zone. In the FeCZ context the situ- 
ation is different, in that only a small fraction of the stellar heat 
flux is carried by convection, and that the temperature gradient 
is already significantly above adiabatic. It is plausible therefore 
that a fossil field of a few kG could indeed suppress convec- 
tion, which would have consequences on observables connected 



to the presence of subsurface convection (e.g. microturbulence, 
C09). One can make an analogy here with the intermediate-mass 
stars (A and late B), where fossil fields do appear to suppress or 
at least significantly reduce the weak helium-ionisation convec- 
tion at the photosphere which manifests itself as microturbulence 
(D. Shulyak, priv. comm.) 

C09 proposed the FeCZ as a possible unifying mechanism 
for the formation (or seed) of different observational phenom- 
ena at the surface of hot, massive stars. The presence of turbu- 
lent subsurface convection appears to be able to trigger at least 
three different classes of physical phenomena: runni ng waves (g- 
modes and p-modes, e.g. lGold reich & Kum arl 19901 C09), solar- 
like oscillations (e.g., C09:lBelkacem et al. 2010) and magnetic 
fields (C09; this work). Each of these perturbations can, in prin- 
ciple, results in observable effects at the surface and/or in the 
stellar wind. Estimating the amplitude, length-scale and time- 
scale of the surface perturbations produced by each mechanism 
is of primary importance to establish the connection between 
FeCZ and observable phenomena. For example wind clumping 
could be produced by small-scale velocity and/or density fluc- 
tuations at the stellar surf ace. These perturbations co uld excite 
the line-driven instability dRunacres & Owo cki 2002) closer to 
the stellar surface than originally predicted from the theory. This 
would reconcile theory w ith observations of the radial stratifi- 
cation of wind clumping dPuls et alj|20 06). However both run- 
ning waves and localized magnetic fields could be responsible 
for such a perturbation at the base of the wind. To understand 
which one dominates, hydrodynamic simulations of the line- 
driven instability are required. These calculations need to in- 
clude, as boundary conditions at the base of the stellar wind, 
the different velocity and density perturbations predicted in the 
case of running waves and surface magnetic fields. While simple 
estimates for these perturbations exist (C09 and this paper), real- 
istic MHD simulations of the FeCZ and the radiative layer above 
it are necessary. Understanding the exact role and properties of 
each perturbation is also fundamental to explain the puzzling co- 
existence of small- and large- scale phenomena observed at the 
surface and in the wind of hot, massive stars (e.g. microturbu- 
lence/macroturbulence and wind clumping/DACs). An effort is 
clearly required to improve existing MH D calculations of sub- 
surface convection dCantiello et al ll2010l) . 

The presence of surface magnetic fields might have impli- 
cations for the evolution of OB stars. This is because, for wind 
confinement parameters r\ <: 1 , the magnetic field can in prin- 
ciple affect the stellar wind mass-loss and angular momentum- 
loss rates. For large scale fields at the s urface, the impact on 
angular-momentum l oss has been studied (Webe r & Davisl 19671 ; 
lud-Doula et al.l 120091) . In the case of the magnetic star <x Ori E 
the theoretical spin-do wn is in agreement with observations 
dTownsend et alj |2010|) . The effect of large scale fields on the 
angular momentum ev olution of main sequ ence massive stars 
has been investigated by Mev net et al.l (1201 ll) . Their results seem 
to show that only fields above a few hundred gauss can have a 
substantial impact. For a given surface field amplitude, higher 
order multipoles (i.e. less coherent fields, as expected in the 
case of dynamo action in the FeCZ) result in a reduced angu- 
lar momentum-loss. Moreover in the case of small-scale fields 
a lot of the flux may reconnect and disappear fairly close to the 
surface, below where most of the wind acceleration takes place. 
Therefore we expect magnetic fields produced in the FeCZ to 
have only a modest effect on the angular momentum evolution 
of massive stars. Nevertheless further investigation of their ge- 
ometry and amplitude is required to confirm this statement. 
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6. Conclusions 

We have shown for the first time that emerging magnetic fields 
produced in the FeCZ of hot, massive stars could reach the stel- 
lar surface. The surface fields expected are localized and with 
amplitudes smaller th an a few hundred gaus s. 

As pointed out bv lHenrichs et al.l d2005l) . a number of unex- 
plained observational phenomena in hot, massive stars are con- 
sistent with surface magnetism. These phenomena include the 
photometric variabi lity and X-ray emission discussed above , line 
profile variability (iFullerton etalJll996t iMorel et al.1l2004 and 
discrete absorption com ponents (DACs, e.g. Cranmer & Owocki 
1996t lKaperetal.lll997h . The magnetic fields emerging from the 
FeCZ could be widespread in OB stars at solar metallicity, and 
play an important role in some or even all these phenomena. The 
metallicity dependency of the FeCZ is such that these effects 
would weaken and eventually disappear with decreasing metal- 
licity (C09). 

Direct detection of magnetic fields emerging from the FeCZ 
has to overcome the problem that, for complex fields, opposite 
line-of-sight magnetic polarities and their respective opposite 
circular polarization signals do cancel when integrated over the 
stellar disc. However if the star is rotating fast enough, it might 
still be possible to detect such fields us ing the Zeeman Doppler 
Imaging technique (ZDI. ISemell lT989). While direct measure- 
ment of FeCZ-generated magnetic fields is challenging, we ar- 
gue that a viable indirect approach is to use high-precision pho- 
tometry - indeed we show that magnetic fields at the stellar sur- 
face of radiative stars should produce bright spots. The finite life- 
time of the magnetic spots, together with rotational modulation, 
could leave a detectable signature in the light curve of hot, mas- 
sive stars. This k ind of variability seems to be present in some 
COROT targets (iDegroote et al.ll2010t IFapics et al.ll201 ll) . and 
might be used to gather preliminary informations on the mag- 
netic field geometry. These information could be used to assess 
the feasibility and the requirements of a targeted direct detection 
campaign. 
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